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Abstract: We study the possibility that dark matter is a baryon of a new strongly 
interacting gauge theory, which was introduced in the low energy theory of Cosmolog- 
ical SUSY Breaking (CSB). This particle can fit the observed dark matter density if 
an appropriate cosmological asymmetry is generated. The same mechanism can also 
explain the dark/baryonic matter ratio in the universe. The mass of the dark matter 
particle is in the multiple TeV range, and could be as high as 20 TeV. 
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1. Introduction 

The most common particle physics models for dark matter involve weakly interact- 
ing particles. They can be broadly classified as WIMPS or axions, with the theoreti- 
cian's favorite WIMP being a neutralino of the Supersymmetric Standard Model (SSM). 
Within string theory, the physics of both of these candidates is closely connected to 
SUSY breaking, because string theory axions generally arise from moduli fields, whose 
mass is related to a superpotential on moduli space. 

One of the authors has recently introduced a new model for SUSY breaking, which 
has no candidate for either WIMP or axion dark matter [1]. The model is based on the 
principle of Cosmological SUSY Breaking (CSB): 

• The (positive) cosmological constant (c.c.) is a discrete tunable parameter, gov- 
erning the number of states in the Hilbert space of quantum gravity in de Sitter 
(dS) space. 

• As the c.c. vanishes, SUSY is restored, with the relation m 3 / 2 ~ A 1 / 4 between 
the gravitino mass and the c.c. A discrete Z n R symmetry is restored in the 
same limit, explaining, in low energy terms, the vanishing of the c.c. in the SUSic 
limiting theory. The limiting theory must have a compact moduli space, in order 
to guarantee that the dS state of the low energy effective field theory is stable. 
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• SUSY breaking is spontaneous in the low energy effective theory, but is induced 
by R breaking terms in the Lagrangian which have no low energy explanation. 
The coefficients in these terms are tuned to guarantee the CSB scaling relation 
between m 3 / 2 and A. 

As a consequence of the first requirement, the low energy effective field theory of 
CSB must contain a goldstino field: a linear supermultiplet which is massless in the 
SUSic, R symmetric limit. In [1] this was taken to be a chiral superfield G, with R 
charge 0. If there are no fields of R charge 2 mod n in the low energy theory, then 
G is naturally massless. R charges were assigned to standard model fields in a way 
that insured the absence of all baryon and lepton number violating dimension 4 and 
5 operators, apart from the term n^H^LiLj (which gives rise to neutrino masses). 
The generation of this term, and of the texture of Yukawa couplings is imagined to 
have to do with physics at the unification scale. There is also an ordinary discrete 
symmetry J 7 , under which G transforms. JF allows the coupling g^GHuHd but forbids 
the conventional /i term. G a is the lowest order T invariant monomial in G. 

High energy physics supplies us with a term MpA 1 ^ 4 f(G/M P ) which violates R 
and implements CSB. The dimensionless coefficients in the function / are tuned to 
guarantee that the c.c. is indeed A. For phenomenological reasons, one must also add 
terms 



We have used an unconventional notation where a lower case label s for a chiral super- 
field S stands for S/M\. The Kahler potential depends, of course, both on chiral fields 
and their conjugates. The functions K and Za are imagined to emerge from integrating 
out degrees of freedom at a scale Mi <C M v <C M P) whose value is determined by RG 
flow in the limiting A = 0, theory. They can be chosen to satisfy all phenomenological 
requirements if M 1 ~ 1 TcV. It is easy to invent strongly coupled theories Q which 
could give rise to all the required properties save one. There is no known example 
of a theory which preserves the R symmetry, and leaves exactly one effective chiral 
superfield which could play the role of G. We will leave this problem to future work 
and concentrate on the problem of dark matter. 

If the coupling functions Za were forced to be logarithms by an accidental U(l) 
with standard model anomalies, then the real part of G could be a QCD axion. How- 
ever, it would have a range of axion couplings ruled out by beam dump experiments. 
Consequently the model has no axion candidates. The basic setup of CSB contradicts 




and, 
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the idea of SUSY neutralino dark matter. The gravitino is the LSP in the CSB sce- 
nario, and its longitudinal components are relatively strongly coupled, so the NLSP is 
not cosmologically stable. 

The only plausible dark matter candidate in this scenario is what we will call a 
Q baryon. That is, we assume the strongly interacting Q sector has an accidental 
symmetry, which renders the lightest particle carrying some accidental U(l) quantum 
number, cosmologically stable. In this paper, we will explore the idea that the dark 
matter is in fact a baryon of a strongly interacting sector with an RG scale of order 
M\. We will see that under a variety of assumptions about the production of this 
particle, this hypothesis is consistent with conventional cosmology. It has the added 
virtue of correlat ing the coincidence between the dynamical scale Mi and the CSB 
scale \J (A 1 / 4 Mp) to the existence of galaxies. That is to say, we imagine that the 
limiting model calculates the value of the scale Mi and the other parameters of e.g. 
the inflaton field, in such a way that the density of Q baryons coincides with what we 
know about dark matter density from observations. Now consider the model of CSB, 
with various values of A. The only values which will produce a model with galaxies will 
be those which satisfy Weinberg's bound. At least within a few orders of magnitude, 
this matches the scale of CSB to M x and the dark energy density to the dark matter 
density (cosmic coincidence). 

We will also see that there is a variety of thermal histories for the universe in 
which Q baryons can be dark matter only if there is a CP violating Q baryon number 
asymmetry. We might imagine a model in which Q and ordinary baryon asymmetries 
were produced by the same mechanism, perhaps explaining the dark/baryonic matter 
ratio of the universe [2] [3]. 

The Hess telescopes [4] have seen a photon signal from the center of the galaxy, 
which might be consistent with a dark matter candidate of mass 15 — 18 TeV, if dark 
matter in the galaxy follows the profile predicted by [7]. It is very hard to find a neu- 
tralino model which can produce such a large mass, basically because weak annihilation 
cross sections decrease with mass. On the other hand, strongly interacting particles 
have mass independent annihilation cross sections and can easily fit this data. 

In the next section, we estimate various cross sections for baryon like objects, 
using large N QCD as a paradigm. The Q theory must differ from QCD since it 
preserves chiral symmetry and is supersymmetric. Nonetheless, we hope that these 
estimates give us a rough guide to the scales involved. We then go on to estimate 
the mass, cross section and primordial asymmetry for which a Q baryon could be dark 
matter. We consider two scenarios: a standard thermal relic abundance calculation, 
and a particular non-thermal production scheme. We find that for reasonable values 
of parameters, the model can fit the data, and perhaps reproduce the Hess signal. To 
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answer the latter question in more detail, one must perform a detailed estimate of 
the photon spectrum one would get from annihilation processes involving a strongly- 
interacting dark matter candidate. We are not sure that the model used by the Hess 
collaboration in order to extract the parameters of a hypothetical dark matter particle 
from their signal, takes into account the physics of a strongly interacting particle. 

We should emphasize that despite our original motivation, our calculations would 
be applicable to any dark matter candidate with new strong interactions of the right 
scale. In particular, we note that our model for dark matter is similar to the hypothesis 
that dark matter is a techni-baryon[5][6]. 

2. Annihilation Cross Sections for Dark Matter With New Strong 
Interactions 

The nucleon anti-nucleon annihilation cross section is usually written in units of the 
pion Compton wavelength, because this is the range of nuclear forces. In fact, this 
parametrization is singular in the chiral limit, when the pion becomes a Goldstone 
boson. It is not correct that the cross section blows up in this limit. 

A better estimate is obtained by thinking about chiral soliton models of the nucleon[8] . 
In such models the nucleon is realized as a classical solution of a large N effective ac- 
tion. The effective Planck constant of this action is of order N, and the scale over which 
solutions vary is the QCD scale. Although these models use the spontaneously broken 
chiral symmetry of QCD in an essential way, they give the same order of magnitude 
results one would expect from general large N considerations. We expect the size of 
a general large N soliton to be given by such an N independent scale, and large N 
soliton masses will be of order N. 

The soliton-anti-soliton annihilation cross section will be given by its classical size 
a ~ Ag 2 and will be more or less energy independent in the regime of interest, because 
the cosmological velocities of these heavy particles will be low. Note that this is not s- 
wave annihilation. The typical orbital angular momentum involved in these collisions is 

of order j^-, where T is the temperature at which the annihilation takes place. Note 
also that the thermally averaged cross section < o~v >, which appears in cosmological 
Boltzmann equations, will be 0{T /iiig) 1 ^ 2 ■ We believe that this is the correct scaling 
even for ordinary baryons, and that conventional calculations of the relic baryon density 
in a baryon symmetric universe are not quite correct. However, this does not change the 
qualitative conclusion of those calculations, namely that we need a baryon asymmetry 
to account for the observed baryon number density of the universe. 
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We note that the reason that we are interested in large N counting is the combina- 
tion of the Hess data, and the constraints on Ag from supersymmetric phenomenology. 
The latter prefers a scale Ag ~ 1 TeV, in order to accommodate the bounds on charged 
superpartner masses, while the former indicates a mass around 15 — 18 TeV for the 
dark matter particle. In a large N model, the baryon mass would be NaAg with a a 
number of order 1 (a ~ 2 in QCD). Thus, we would want Na ~ 15 — 18. These are not 
unreasonable values. For example, the best of the inadequate models for the Q theory, 
studied in [1] was an SU(4) SUSY gauge theory. For N = 4, we require a ~ 4, about 
twice the value in QCD. 

We emphasize however that we do not know the details of the model which the Hess 
collaboration used in quoting 15 — 18 TeV for their best fit to the dark matter candidate. 
In particular, for weakly coupled neutral dark matter, the direct photon annihilation 
signal is suppressed by a power of (a/V) 2 relative to photons produced from decays 
of particles with direct coupling to the dark matter. There is no such suppression for 
strongly interacting neutral composites of charged particles. For example the large N 
nucleon magnetic moment is order e(= y/Ana^jN in Aqqd units. Hess has not yet 
seen the characteristic turnover in their photon signal, which would be expected from 
dark matter annihilation, and the question of astrophysical explanations for the signal 
from the galactic center is still controversial. It is perhaps premature to try to fit their 
spectrum. 

However, it is clear that in order to really confront an eventual dark matter signal 
from Hess data, we need a much better estimate of the photon spectrum produced by 
a Q baryon. In addition, since we find that for most values of the reheat temperature 
of the universe, we must invoke a Q baryon asymmetry to account for the observed 
dark matter density, the annihilation signal will be proportional to the small density 
of anti-£ baryons. We have not yet done the calculations to determine the range of 
parameters for which we would expect a significant annihilation signal from the center 
of the galaxy. In the rest of this paper, we will choose an annihilation cross section of 
order Ag 2 and parametrize our results in terms of the Q baryon mass mg > Ag, Ag, 
and an asymmetry. 

Our description of the Q baryon will utilize the following characteristics of a soliton 
model: energy independent annihilation cross section much larger than the scale of its 
Compton wavelength, and thermal production at energies well below its mass. The 
latter is a well known characteristic [9] of solitons in weakly coupled field theory. Finally, 
we will parametrize the Q baryon mass as NaAg, with Ag ~ 1 TeV, in order to suggest 
the large N scaling of soliton masses in strongly coupled gauge theories with large 
gauge groups. 
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3. The Relic Abundance of Q baryons 



We will denote by VIq the fraction of the observed density of the universe in Q baryons 
plus anti-baryons. To match the observed dark matter abundance, we require Qg = 
— = .24, using the data from WMAP which specifies Vt m = .29 ± .07 and fif, = 

Per 

.047 ± .006.(10] If ric is the number of Q baryons per comoving volume, then this can 

be written Q — "? mg — n ° ms 

De writxen \ i - 3H 2 /8wG - im4h 2. 10i ^v ■ 

J cm* 5 

Writing today's value of the Q baryon abundance (the ratio of the number of Q 

?y) Yo = 

s Y m g 



baryons per comoving volume and the entropy) Y = this condition becomes 



n = .24 



1.054/i 2 • 10 4 ^C 



Thus we require Y = ^ . We will write mg = Na TeV, treating 1 TeV as the analog 
of the QCD scale for the Q gauge theory, and applying a large N scaling rule for baryon 
masses. In QCD N = 3 and a ~ 2. Our point is that the analog of a baryon mass 
could be quite a bit higher than 1 TeV. For example N = 6 and a ~ 3 would give us an 
18 TeV dark matter candidate, as would be required by the interpretation of Hess data 
in terms of dark matter annihilation. With this parametrization, the required value of 

A in- VJ 

the abundance is Y ~- ;:! 



Na ■ 

The relic abundance of Q baryons depends on some assumptions about the evolution 
of the universe at the TeV scale and above. We assume that there was a reheating 
process which gives rise to a radiation dominated universe at some temperature Tr#. 
This might be due to primordial inflaton decay, or the later decay of some other massive 
particle which dominates the energy density before it decays. We call the width of the 
particle T x ■ If T RH > 1 TeV, the Q gauge theory is thermalized by X-decay and the 
post-decay distribution of Q baryons is given by the thermal ensemble. Note that this 
is true even when mg T RH . In this regime of parameters, the Q baryon is a thermal 
relic, and we find that, in the absence of an asymmetry, the relic abundance is too 
small to explain the observed dark matter density. 

For T RH < 1 TeV, Q baryons are produced non-thermally and we must be a bit 
more specific about the dynamics. For a weakly coupled X particle, mx 3> T R h and 
we can still have X decays into Q baryons. Suppose first that mx 3> vrtg so that we 
can treat the Q baryons as just another massless species. If we assume the couplings to 
Q baryons are not suppressed relative to standard model particles we get a branching 
ratio of order 10~ 2 into Q baryons. The decay will be reasonably rapid, so we neglect 
annihilation processes during the decay period and obtain an initial abundance of 

Yo ~ lO" 2 ^. 

mg 
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If the Q baryon were massless, this ratio would just be the branching ratio 1CT 2 . The 
additional suppression is our estimate of the number of Q baryons per photon that 
result from the thermalization process. 

Throughout the interesting range of parameters, the X particle life-time is short 
enough to neglect annihilation in the calculation above. Now we can evolve the result- 
ing Q baryon densities according a Boltzmann equation driven only by the annihilation 
of g and g. As in the discussion of solitonic dark matter abundance in Griest and 
Kamionkowski [11], the thermally averaged annihilation cross sections have tempera- 
ture dependence given by : < a\v\ >= °"o("^ L ) 1 ^ 2 - Also, we will assume that T RH < mg. 
In this case there can be no process in the Boltzmann equation that creates Q baryons 
because it is not energetically favorable. The Boltzmann equation for the evolution of 
Q baryons is: 

h g + 3Hn g = — < a\v\ > n 2 

Letting Y = n g /s and x = mg/T we get 

dY = x 1 / 2 a sm p i y2 
dx 1.67gt /2 m 2 g 

Since s = ^g* s T 3 , we can then write: 

dY kY 2 



dx x 5 / 2 



where k = mg2n a ^ m ^ . 

1.67s;' 45 



Here we will assume an average g* m g* s pa 50. 
Defining an order one parameter (5 such that o"o = jjyf^yjz, k ~ 4.5 • 10 15 Na/ [3 2 
The solution to this equation is: 



Y, 



1 



final i 2k ( 1 



J_ -I- 2fc / 1 1 \ 

Yi ^ 3 ^1372 -J72) 



x f 



Notice a few properties of this solution. The present day temperature is so low 
that ~ 0. Hence either the ^ term or the term dominates, depending on 

Xf r * x i 

T RH . The y: term dominates for T RH < .3Na/(3 2 MeV. A reheat temperature in this 
range would be inconsistent with nucleosynthesis, so we can ignore this term. Thus, Yf 
is determined by: 

f 2k 
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where x { = (m g /T RH ). 

In a general model where we do not fix the mass of the Q baryon or the exact cross 
section, we can get an upper bound on Trh from our requirement that Y = AA '^ a — : 

Trh > .WSp^NaTeV 1 

For reheat temperatures below this value, the Q baryons will dominate the universe. 
Thus we find a small window 1 > > .008/3 4 / 3 iVa, where non-thermal, symmetric 
Q baryon production could account for the observed properties of dark matter. In 
particular, for typical values Na ~ 10 and (5 ~ 1, we find that this window has a width 
of about an order of magnitude. However, this range for the reheat temperature does 
not conform to our prejudice that mg is substantially larger than Ag. We also note that 
there was no loss of generality in our assumption that m x » mg. If this assumption 
is not valid, then Trh is quite low, and Q baryons would be overproduced as long as 
m x > mg. 

For Trh > 1 TeV, the thermal relic abundance is too small to account for the 
observed dark matter, but we can remedy this by postulating an asymmetry. The 
simplest possibility is that the asymmetry is generated directly in the decay of the X 
particle, in which case we have the standard result that 

Y = e G , 

m x 

where 

v x (x -^f)-r x (x — j) 

f 



T is a decay rate, / and / are all possible final states, and Bj is the total Q baryon 
number of the final state /. is the corresponding asymmetry in ordinary baryon 
number. In order to match the observed dark matter density and the observed baryon 
density, we need 

^ « J- x 10- 2 , 
e B 2Na 

and 

T 

e B — w 8.6 x 10 -11 . 

m x 



1 When this number is larger than 1 TeV the calculation is not self consistent, because Q baryon 
production is thermal. 
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In our model ~ v/m j ^ mp ) 2 is bounded from below by the requirements that 

the X couplings to ordinary matter are at most Planck suppressed, and that the X is 
massive enough to produce the Q baryon in its decays. Thus ^ K > (^p) 1 / 2 x 1CT 3 . 

We see that the e parameters must be very small in order to account for the observed 
asymmetries. In fact, small baryon number violating branching ratios arise naturally 
if we assume that X is a slow roll inflaton, /, with a "natural" potential of the form 
fi 4 f(I/mp). A theorem of Nanopoulos and Weinberg[12] tells us that asymmetries can 
arise only at second order in baryon violating couplings. Let us assume the decay of the 
inflaton is mediated primarily via dimension 5 operators. Then even if the dimension 
five couplings involve CP violation and baryon number violation, we will find that 
e_B ~ (— ) 2 . ^y e a j go faoyQ the order of magnitude estimate ~ (ullU/ 2 so ^hat 

Yb ~ (^) 5/2 . 
m P 

This will fit the observed baryon asymmetry if 

mi ~ 10~ 4 mp 

Note that this gives an inflation scale /j, close to the unification scale. 

In this context we might attempt to explain the further suppression ^ ~ 1CT 3 
by postulating that (perhaps as a consequence of the R symmetry introduced in [1]) 
the leading contribution to the ^-baryon asymmetry comes from the interference of a 
dimension 5 and dimension 6 coupling of the inflaton, and is suppressed by a further 
power of — . This is off by a factor of 10 but our estimates are so crude that we can 

r mp J 

consider this a success. 

Indeed, we proposed this simple model not because we think it has to be right, but 
to show that reasonable calculations of both the dark matter and baryon abundances 
can be obtained for our new form of dark matter. 

To summarize, we probably need asymmetric production of Q baryons to make them 
an acceptable dark matter candidate. We outlined a plausible model of asymmetric 
production in inflaton decay, which could naturally explain both the baryon asymmetry 
of the universe and the dark matter density. 



4. Conclusions 

We have shown that a baryon-like state of the new, strongly interacting, Q theory, 
which was introduced in [1] to implement Cosmological SUSY breaking, is a promising 

2 M is the scale of irrelevant couplings of the X particle to the standard model, which are responsible 
for the decay. 
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dark matter candidate. The new strong interaction scale is around 1 TeV and the Q 
baryon mass is somewhat higher, perhaps as high as the 15 — 18 TeV needed to fit the 
Hess data on photons from the center of the galaxy, if the explanation for that data 
turns out to be dark matter annihilation. We saw that this sort of baryon to interaction 
scale ratio was natural in the context of large N scaling with N ~ 5 — 7. While the 
Q theory is probably not as simple as SU(N) QCD, there is reason to believe that 
reasonably large baryon masses are a more general phenomenon. 

This sort of dark matter candidate allows one to contemplate a simple explanation 
of the dark matter to baryon ratio, since the asymmetries in baryon and Q baryon 
number might have the same physical origin. We need to explain a factor of order 
10 3 in these asymmetries, in order to fit the data. We constructed a plausible model 
in which both asymmetries are generated in inflaton decay. To explain the size of 
the asymmetries we invoked the R symmetry of [1] and dimensional analysis. There 
are more scenarios for baryogenesis in the literature than there are authors on this 
paper, and it is entirely plausible to us that a more elegant mechanism could be found. 
However, our simple model might work, and it might be the right answer. 

Much more work needs to be done to sort out signatures of such hyper-strongly 
interacting dark matter, as well as to explore a variety of models for the production of 
baryon and Q baryon asymmetries. In addition, it will be necessary to find out more 
about the dynamics of the as yet mysterious Q theory, which gives rise to these new 
particles. 
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